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GPS RECEIVER TRACKING SYSTEM 
CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims the benefit of U.S. provisional application No. 

60/256,870, filed on December 20, 2000, which is hereby incorporated by reference in its 
entirety. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0002] The present invention relates to the field of Global Positioning System (GPS) 

receivers receiving transmission GPS signals from satellites and processing the GPS signals 
for obtaining information, and storing the information. The present invention particularly 
relates to a GPS receiver tracking system tracking a high dynamic vehicle in accordance with 
GPS data received from a translator placed on the high dynamic vehicle. The system is 
preferably implemented entirely using software. 

2. Description of the Related Art 

[0003] The Global Positioning System (GPS) is a sateUite-baised navigation system 

fliat continuously transmits timing, frequency and satellite position information to potential 
users. The GPS consists of a full constellation of four satellites in half geo-synchronous 
orbits. The GPS satellites continuously emit coded GPS signals The GPS signals can be 
received Scorn at least four GPS satellites at any point on or near the Earth. 
[0004] A GPS signal contains timing information that allows a user to detemiine the 

time elapsed for the GPS signal to transverse the distance between the GPS satellite and a 
receiver receiving the signal. By knowing the time the GPS signal left the GPS satellite, fiie 
time the GPS signal arrived at the user, and the speed of flie GPS signal, flie receiver can 
determine the distance from itself to the GPS satellite. The measured range is referred to as 
**pseudorange" because there is generally a time difference or offeet between timing clocks on 
the satellites and the GPS receiver clock. Thus, for three dimensional position 
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determinations, at least four sateUite signals are needed to solve for four unknowns, i.e., the 
time offiet and the three dimensional positions of the satellites. By knowing the orbital 
position of four GPS satellites (ephemeris data), and the distance from itself to each of four 
GPS satellites, the receiver can successfully triangulate its own position. Receiver position 
calculations performed using greater than four satellites generally have an improved accuracy. 
[0005] The GPS signal emitted by the satellites contains an L-band carrier component 

(LI) transmitted at a frequency of 1575.42 MHz. The LI carrier component is modulated by a 
coarse acquisition (CVA) pseudorandom (PKN) code component unique to the satellite .anda 
data component. The PRN code provides timing information for determining when the GPS 
signal was broadcast and identifymg which satelUte emitted the signal. The data component 
provides information defining the satellite q)hemeris, satellite clock corrections and other 
GPS information. The carrier component aUows a receiver to more easily acquire the GPS 
signal. 

[00061 It is known in the art to use Ihe Global Positioning System (GPS) for testing 

and evaluating purposes, such as for testing and evaluating the performance of moving 
vehicles. However, accurate development of test data for hi^ dynamic vehicles has proved 
problematic, as well as for test data obtained from multipath GPS signals and GPS signals 
reflected from the ocean. 

[0007] GPS receivers have beai developed for testing of high dynamic vehicles. 

However, the receivers are limited for use with very high dynamic vehicles, such as missfles 
having a short time of flight and small dimensions. GPS receivers are generally too large to 
be fitted on a small missile, and flie time needed for obtainmg a first set of data may exceed 
the «itire flight time. 

[0008] Typically, use of GPS receivers for platforms such as small missiles has been 

replaced by Ihe use of translators. A translator is a device capable of downlinking the 
complete GPS spectrum to ground recording equipment for post flight tracking and trajectory 
reconstruction. The translators have applications for higji dynamic guidance, range safety and 
post-test trajectory analysis. However, translators currently available tend to be costly and 
highly proprietary, limiting the possibility of modifying the translator for special appUcations. 
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Another drawback to available translators is the need for accessing a large bandwidth that is 
not always available at test ranges. 

[0009] In one suggested solution, recorded data obtained from the translator and 

recorded data obtained by a groimd receiver approximately simultaneously with the data 
obtained by the translator are processed by filtering and smoothing algorithms to produce a 
smooth trajectory. However, during the boost phase of the fligjit, the signal-to-noise ratio 
(SNR) is very low. 

[0010] Accordingly, there exists a need for a method and a system for processing 

recorded GPS data obtained form a translator mounted on a mobile dynamic vehicle and 
recorded data obtained approximately simultaneously witii the data obtained by the translator 
by a ground receiver, in which the tracking results are accurate for all phases of movement of 
the vehicle, where the vehicle is a high dynamic vehicle. 

SUMMARY OF THE INVENTION . 

[0011] It is an object of tiie present invention to provide a method and a system for 

processing recorded GPS data obtained fiom a translator and recorded data obtained 
approximately simultaneously with the data obtained by the translator by a ground receiver, in 
which the tracking results are accurate for all phases of the flight of a very high dynamic 
vehicle or for GPS signals reflected off the water. 

[0012] It is a fiuther object of the present invention to minimize processing time while 

maintaining accuracy wheii processing recorded data obtained from a translator and recorded 
data obtained approximately simultaneously with the data obtained by the translator by a 
ground receiver. 

[0013] Yet a further object of the present invention to provide a GPS receiver which 

performs tracking and recovery of timing information using software instead of a combination 
of hardware and software. 

[0014] In keeping with these and ofhCT objects of the present invention, a receiving 

station tracking system for receiving GPS data and tracking a mobile vehicle is provided. The 
receiving station tracking system includes at least one processing stage receiving and 
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processing at least a portion of the GPS data, and a control module receiving and filtering the 
processed GPS data and outputting at least one control signal in accordance with filtering 
results. The at least one control signal controls decreasing or increasing a quantity of 
processing stages of fee at least one processing stage, 

PtT?TRP DESCT^TPTTON OP TWK TO? AWINGS 

[0014] FIG. 1 is a schematic diagram of a prior art recaving station trackmg a missile; 

[0015] HG. 2 is a block diagram of an exemplary receiving station in accordance with 

the present invention; 

[00161 HG. 3 isablockdiagramofatrackingmoduleofthereceivingstationin 

accordance with the presait invention; 

[00171 FIG- 4 is a block diagram of a translator .tracking module of the tracking 

module in accordance witti the present invention; 

[0018] HG. 5 is ablock diagram of an exemplary processing stage of the translator 

tracking module in accordance with tiie present invention; 

[0019] FIG. 6A is a plot of experimental results for an in-phase signal and noise 

versus time; 

[00201 PIG. 6B is a plot of experimental results for a quadrature signal and noise 

versus time; and 

[00211 FIG. 6C is a plot of experimental results for a spectrum of signal and noise 

versus firequmcy. 

DFTAn.ED DESCRIPTION OP THE PRP FFT^T^ED EMBODIMENT 
[0022] In FIG. 1. a missile 10 equipped with a conventional translator 14 is shown in 

flight while being tracked using GPS. As is known in ihe art, during the flight of the missile 
10, the translator 14 receives a GPS signal including an LI carrier firom at least four sateUite 
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vehicles (SVs) 18n, where n=l-24, of 24 available SVs 18n. The translator converts each 
received signal to an S band signal (approximately 2265.5 MHz), and transmits the S-band 
signal to a nearby ground receiving station 22. The receiving station 22 receives the S band 
signal, converts it back to a GPS signal and processes the GPS signal and records the 
processed GPS signal as translator GPS data. 

[00231 During the flight of fte missile 10, the receiving station 22 receives a GPS 

signal including an LI carrier from at least four SVs 18(n) (the same set of SVs 18(n) due to 
the proximity of the missile 10 during flight and the receiving station 22). The receiving 
station 22 and the translator 14 receive the GPS signals simultaneously or approximately 
simultaneously. The receiving station 22 processes and records the GPS signal as receiver 
station GPS data. The receiving station 22 performs processing on the recorded translator 
GPS data and the receiving station GPS data. 

[0024] During the processing, described in furflier detail below, a saies of Fast 

Fourier Transforms (FFTs) are repeatedly performed on a set of sampled translator GPS data 
for correlating the translator GPS data. After each repetition the results of the FFT processing 
and the degree of correlation are analyzed and used to reduce processing elements by 
reducing the set of sampled translator GPS data for the next repetition. Reduction of the 
number of processing elements allows the algorithm to run faster. 

[0025] The object being tracked may be any vehicle located proximate flie receiving 

station 22 and equipped with the translator 14. In the example provided in FIG. 1, flie object 
being tracked is a missile 10, such as a Stinger, Hydra rocket or Tomahawk, having high 
dynamic flight qualities in which there is a high variation of velocity and acceleration values. 
During the boost stage of flight, the acceleration of the missile begins at nil, rapidly climbs to 
a high peak value, and continues to rapidly change. During fte coast stage of fli^t the 
acceleration ofthe missile changes gradually. During the landing stage of flight fee 
acceleration drops to nil. 

[0026] Translators are known in fee art. The translator 14 preferably multiplexes 

inputs ftom up to four dual frequency antennas, downconverts feem to an intermediate 
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frequency (IF) where they are filtered and cjombined with a pilot cairier. The signal is then 
upconverted to S-band and transmitted. 

[0027] FIG. 2 shows a block diagram of an exemplary receiving station 22. The 

receiving station 22 includes at least one antenna assembly (not shown) for receiving signals 
transmitted by SVs 1 8 (n) and flie S-band signal transmitted by the translator, and a 
conventional S-band receiver 210. which receives the S-band signal transmitted by the 
translator and converts it into a GPS signal. The receiving station 22 further includes a GPS 
downconverter, sampler and formatter (DSP) 214 that receives the GPS signal converted from 
the received translator S-band signal and the GPS signal received from the SVs 1 8(n). 
Preferably, the GPS DSF 214 samples the receiving station GPS data and translator GPS 
data at the same rate. 

[0028] The GPS DSF 214 downconverts the received receiving station and translator 

GPS signals each into a pair of in-phase (I) and quadrature (Q) signals, samples the signals at 
approximately 23.333 MHz, formats them into 32-bit words, smd passes them through a 32-bit 
high speed I/O 218 to a PCI bus 222 as tune ordered signal pairs I(tO, Q(ti), where tH(At), 
and At=a tune interval per input data sample. A computer 226, i.e. a computer capable of 
running a network operating system, such as Windows NT ™, controls the PCI bus 222 
activity and routes the sampled data to a commercially available higb-rate data recording 
system 230 having a software interface. The recording system 230 places the sampled data on 
a hard drive farm 234 in a proprietary format. The computer 226 may access the sampled data 
using the software interface of the recording system 230, format the sampled data into 
standard network operating system files, and store them on an array of removable SCSI 
peripheral hard drives 238 accessible to the computer 226 via a SCSI adapter 242. 
[0029] The receiving station 22 fiirlher includes a tracking module 244 for processing 

the sampled data associated wilh the receiving station GPS data and the translator GPS data. 
The trackmg module 244 accesses the sampled data stored on the hard drives 238. The 
tracking module 244 may access the sampled data via the computer 226 as shown, or the 
tracking module may include necessary software and hardware for accessing the hard drives 
238. The tracking module 244 may fiirfher be included within the computer 226 or a dififerent 
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computer. It is to be understood that the receiving station 22 of FIG. 2 is exemplary, and that 
the receiving station GPS data and the translator GPS data may be processed and stored by 
other means so that the data is stored as time ordered sampled digital data and is accessible by 
the tracking module 244. 

[0030] FIG. 3 shows the tracking module 244. The modules of the tracking module 

244 may be implemented by software, hardware or a combination thereof. The tracking 
module 244 includes a user interfece UI 3 10, a receiving station data tracking module 314, 
aides data database 3 1 8, and a translator data trackmg module 322. The UI 3 10 includes 
means for exchanging data wife a xiser, such as via a keyboard, joystick, microphone, touch 
screen, etc., where the user can input information, and a GUI, speaker, etc., where the tracking 
module 244 can ou^ut information to the user. 

[0031] The receiving station data tracking module 3 14 uses tracking and/or Doppler 

aides, performs a sky search wifeout using aides, uses other known methods and algorithms 
for processing the receiving station GPS tune ordered sampled data or a combination thereof, 
in order to extract data from fee receiving station GPS data, such as fee navigation message 
and pseudorange. The tracking and Doppler aides are initial estimates specifying which 
satellite is in view, satellite orbit data and what carrier Doppl^ frequency is expected. The 
tracking and Doppler aides may be provided by user input, remote source input, or accessed 
from fee aides data database 318. The aides data database 318 stores known information about 
fee GPS signals transmitted by fee SVs 18(n): Once data, such as fee navigation message and 
pseudorange values data, is extracted from fee receivmg station GPS sampled data, selected 
data such as navigation message bits, fee absolute sample index for each C/A code epoch, and 
time of week (TOW) are captured into a capture file to be provided as tracking aides for 
processing fee high dynamic signal from fee translator. ' 

[0032] The translator data tracking module 322 processes fee translator GPS data by 

accessing fee capture file provided by fee receiving station data tracking module 3 14 for 
obtaining pseudorange and Doppler data associated wife fee missile being tracked. 
[0033] FIG. 4 shows fee modules and data flow of fee translator tracking module 322. 

The modules of fee translator tracking module 322 may be implemented in software, 
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hardware, or a combination thereof. The translator tracking module 322 includes a series of 
processing stages 410, a coordinate aides module 414, a local carrier generator 418, a PN 
Code Generator module 422, a filter and search module 426, a stage control module 430, an 
input control module 434, and a peak interpolation module 438. Each processing stage 410 
includes a complex multiply module 450, a delay module 454, a pair of multiplier modules 
458, a pair of sum and dump modules 462, an N-Point FFT module 466, and a magnitude 
module 470. A signal pair I, Q, where I,Q are periodic signals such as sinusoidal waves, for 
each value of a sample time tj, and flie capture ffle are provided as inputs to the translator 
tracking module 314. The translator tracking module 3 14 ou^uts a pseudorange and Doppler 

value for each signal pair I, Q. 

[0034] M processing stages 410(s), where s=l to m are shown in FIG. 4 where a pak 

of translator GPS data input sample signals I,Q for a given time, ti are processed relative to a 
FN code generated by the PN Code Generator 422 and a local carrier signal including periodic 
complex sinusoidal waves (I*(ti), Q*(ti)), for the time ti, provided by the local carrier 
generator 41 8. The PN code is delayed by a different delay at each stage as it passes throu^ 
delay modules 454. Each delay module 454 generates a delay d, where d is typically half a 
chip. In the exaraple shown, at processing stage 410(1) the PN code is delayed by d=0. At the 
processing stage 410(m) the PN code is delayed by d=m-l. As described in furflier detail 
with refwence to FIG. 5, each processing stage 410(s) includes the N-point FFT module 466 
for performing Fast Fourier Transforms on the signals being processed and outputs N 
magnitude values, where N is the number of data points per FFT conversion performed by the 
N-point FFT module 466. 

[00351 The capture file is used for control of Ihe PN Code generator module 422. The 

capture file, or a portion thereof is accessed by the coordinate aides module 414. The 
coordinate aides module 414 receives a signal firom the input control module 434 indicating 
which sample of the series of input GPS samples is currently being hspnt, and coordinates the 
data in die capture file with the input translator GPS sampled data. The coordinate aides 
module 414 fiirthea: recovers a C/A Code epoch index tcom data in the capture file for setting 
the code position (timing) for processing by a reference processing stage 410(s), such as the 
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first processing stage 410(1). The coordinate aides module 414 provides navigation message 
data, including polarity of the PN Code, to flie FN Code generator module 422, and further 
provides a receiving station PN Code Epoch position, which is timing information encoded in 
data of the capture file, to the Peak Algorithm module 438 and to the filter and search module 
426. 

[0036] The local carrier generator 418 provides the local carrier signal (I*(ti), Q*(ti)) 

to the pair of multiplier modules 458, which downconverts the translator GPS signal to near 
baseband The local carrier generator 418 includes is an oscillator, such as a number 
controlled oscillator (NCO). The NCO includes a phase accumulator to which a delta, which 
sets the oscillator firequency, is added for each input sample processed. The output of the 
phase accumulator is used to generate a complex sinusoid used as flie local carrier signal. The 
rate of the NCO of the local carrier generator 418 is optionally controlled by a rate control 
signal provided by the filter and search module 426, or the rate of the local carrier may be 
fixed. Due to a Doppler effect when Fast Fourier Transform processing is performed within 
the processing stages 410, control of flie rate of the local carrier generator 418 is typically not 
needed, but may be adjusted by the local carrier generator 418 when tiie filter and search 
module 426 determines that correlation of input translator data and receiving station GPS data 
is extremely poor. 

[0037J The PN Code generator 422 generates a PN code having a pattrai of bits, 

where the bits form a random number. The PN code generator 422 may include an oscillator 
such as a number-controlled oscillator (NCO). A predetermined PN code is associated with 
each satellite. The navigation message output by the coordinate aides module 414 controls the 
polarity of the PN Code Generator 422, effectively wiping off the navigation message fi-om the 
receiver GPS signal, allowing the sum and dump module pair 462 to accumulate across 
message bit boimdaries without degradation of performance. A phase and rate control signal 
generated by the filter and search module 426 controls the phase and rate of the NCO of the 
PN code generator 422. The phase and rate control signal provided by the filter and search 
module 426 controls the PN code by shifting the PN code and the rate of the code by 
controlling the rate at which the bits fomiing the PN code are transmitted via the NCO. 
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[0038] The filtra- and search module 426 uses tibe C/A code to control setting the C/A 

code position (timmg) for the reference processing stage, i.e., first processing stage 410(1). 
The filter and search module 426 further filters the N magnitude outputs of the M processing 
stages 410(s), so that the filter and search module 426 filters a matrix of N*M values. The 
filter and search module 426 then searches the filtered values for a peak, processes the 
location of the peak and a Peak-to-Noise ratio, and determines if the results of &e processing 
are satisfactory, i.e., if the correlation between the receiving station GPS data and the current 
translator GPS data input sample has been achieved. 

[0039] If flie filter and search module 426 determines that correlation has not yet been 

achieved, the filter and search module 426 then generates control signals based on the location 
of the peak found during the search with reference to the aides in ihe capture file of the 
receiving station, the Peak-to-Noise ratio, and the PN code epoch position signal received 
fiom Ihe coordinate aides module 414. The control signals generated by the filter and search 
module 426 include at least one of the rate control signal to the local carrier generator module 
418, the phase and rate control signal to the PN code generator module 422, a control signal to 
the stage controller module 430, and a control signal to the input control module 434. Upon 
determination of satisfactory processing results, the output of the processing stages 410 is 
provided to the peak interpolation module 438. 

[00401 The filt«- and search module 426 effectively controls the size of the M*N 

matrix of values that it processes based on the results its search for a peak. The matrix niay be 
enlarged by instructing the stage control module 430 to add selected processing stages 410(s), 
or the matrix may be reduced for reducing processing steps. The processing stages 410(s) to 
be added or removed from the series of processing stages 410 are selected in order to position 
the expected peak towards the center of the matrix. Furthramore, the expected peak is seeked 
out by adjusting the phase and rate of the PN Code Generator 422, and the rate of the local 
carrier generator 418. 

[0041] The stage controller module 430 receives the control signal from the filter and 

search module 426 with instructions to add or remove processing stages 410(s). The 
instructions may include which processing stages 410(s) to add or remove. For example, the 
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instructions may instruct the stage controller module 430 to add or remove only the 
processing stages 410(s) where "s" is a selected number of lowest possible values, highest 
possible values, evenly distributed among low and high values, or unevenly distributed values 
of "s". Thus, the filter and search module 426 controls the number and distribution of 
processing stages 410(s) used during processing of each translator GPS data input sample. As 
processing of the translator GPS input data proceeds, the quality of peak detection typically 
improves, and the filter and search module 426 preferably instructs the stage controller 
module 430 to reduce the value of **m", thus reducing the amoimt of processing, and therefore 
the processing time. 

[0042] The input control module 434 receives the control signal from the filter and 

search module 426, accesses the translator station GPS input data samples and sequences to 
the next translator GPS data input sample I,Q for time t+At, where At is the time interval per 
input data sample, in accordance with the control signal. The local carrier generator 418 and 
the PN code generator 422 also advance to the next output sample based on the programmed 
rate. Thus, the filter and search module 426 controls when to proceed to the next translator 
GPS data input samples. Fxxrfhennore, as described above, the input control module 434 
provides the coordinate aides module 414 with a signal indicating which sample of the series 
of input translator GPS data samples is currently being input to the processing stages 41 0(s). 
[00431 The peak interpolation module 438 is used to determine the peak value and 

Doppler value mote precisely. The peak interpolation module 438 processes the receiving 
station PN Code Epoch position and the time of week aides together with the data output by 
the processing stages 410 and filtered by flie filter and search module 426, and outputs the 
pseudorange and Doppler values for the current translator GPS data input sample being 
processed. 

[0044] FIG. 5 shows an exemplary processing stage 410(p) in detail, where p is the 

value of s for a prompt sample, as is known in the art Table 5, below, defines the terms 
shown in FIG. 5. 
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TABLES 

FFT Acquisition of GPS Signals for Processing Stage 410(p) 

= Received and downconverted in phase data samples 
Q(^^) = Received and downconverted quadrature data samples 

u - mi) 

A^; = Time interval per input data sample 

l*(t^ = In phase output of the number - controlled oscillator (NCO) of the local 

carrier generator 418 
Q*(^^ =: Quadrature output of the NCO of the local carrier generator 418 
C(^^ = Cosine output of the phase detector (complex multiply module 450) 

= lifO I*(/e^+Q(^/)Q*W 
S(f^ = Sine ou^ut of the phase detector (complex multiply module 450) 

= l{tO Q*aiJ-Q(^i)i*W .... 
lApif O = Prompt (i.e., on time) output of the local code generator for delay = p- 1 
x^(/y) = Complex FFT input samples, prompt code alignment 

tj = j{At2) 

/M2 = Time interval between FFT input samples = p(/^i) 
^ = Sample rate of the FFT input samples = 1_ 

;;fy = Nyquist frequency = IJs 

2 

irfsffd Span of frequencies that can be unambiguously identified by the FFT 

algorithm 

Jsf = Number of data points per FFT conversion 

T = Time span of data points going into each FFT conversion - N(At2) 

Xp{f0 = Discrete Fourier transform of the prompt input samples 
Ppifki Power of the prompt input signal in the Wh frequency bin 

= {Real[X^(/i)]}'+ {Imag[^Oi)]}' 

M = Width of the FFT frequency bins = i 

T 
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[0045] The complex multiply module 450, the multiplier module pair 458, the sum 

and dump module pair 462, and the square module 470 of the exemplary processing stage 
410(p) operate as is known in the art for wiping off the navigation message from the translator 
GPS data input sample. Translator GPS data input samples, I(ti), Q(ti) signal pair and local 
carrier I*(ti), Q*(ti) signal pair from the local carrier generator 418, respectively, for time t= ti, 
are received by the complex multiply module 450, (also known as a phase detector) of the 
processing stage 410(p). The complex multiply module 450 multiplies the I(ti), Q(ti) and the 
I*(ti), Q*(ti) signal pairs and outputs C(ti) and S(ti) as is known in the art. The complex 
multiply module 450 performs a last stage of downconversion to ordered sequence (I(ti), 
Q(ti)) using the ordered sequence I*(ti), Q*(ti) received from the local carrier generator 418. 
[0046] The C(ti) and S(ti) signals are each multiplied via a respective multiplier 

module of multiplier module pair 458 by the signal Mp(ti), as is known in the art The signal 
Mp(ti) is the PN code output by the PN code generator 422 and delayed by an appropriate 
delay d generated by delay modules 454, where the delay d== p-1 . A first and second sum and 
dump (also known as accumulate and dump) module of smn and dump module pair 462 
respectively receive the output of one of each of the multiplier modules of the multiplier 
module pair 458. The respective sum and dump module pair 458 simis the output of the 
multiplier pair 458 for P points, where P is a predetennined number, and outputs the result x 
R(tj), and xi(tj), respectively, as is known in the art. Timing of the sum and dump module pair 
462 may be independent of the timing of the bits in message bits recovered by the coordinate 
aides module 414. Each input sample (I(ti), Q(ti)) is processed by the complex multiply block 
450, the multiplier module pair 458 and the sum and dump module pair 462. Once an input 
sample is processed, the next input sample is sequenced to, and is processed. When the sum 
and dump module pair 462 accumulates a predetermined number (P) of samples, the output of 
the sum and dump module pair 462 is transferred to N-Point FFT module 466. 
[0047] The local carrier generator 418 together with tiie complex multiply module 

450, multiplier module pair 458 and the sum and dunip module pair 462 efifectively perform 
matched filter processing of the input translator GPS data, as is known in tie ait 
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[0048] The N-Point FFT module 466 receives XRCtj) and xi(tj), and upon accumulating 

N samples from the sum and dump module pair 462 the N-Point FFT module 466 performs a 
discrete Fourier transform of the prompt input samples over N points. The N-Point FFT 
module outputs N complex numbers Xp(fk)=XR(fk)+i*Xi(fk) to respective first and second 
square modules 472 of flie magnitude module 470. The first and second square modules 472 
respectively output {Real[Xp(fiO]}^ and {ImagP^p(fk)]}^ The outputs of the first and second 
square modules 472 are summed by sum module 474 and output as Pp(fk), having N values, 
including one value for each frequency bin of the FFT, effectively perfomiing a power 
spectrum analysis. 

[0049] An example is now discussed of the opemtion of the translator tracking module 

322. m processing stages 410(s) are enabled for processing. The number m and the selection 
of stages are selected either by an algorithm or entered by an operator. Selection of m is 
typically based on known criteria, such as the location of tiie receiving station 22 and flie 
range of the missile 10 or other trajectile. Preferably, m is initially set to allow for a delay 
corresponding to twice the maximum distance between the receiving station 22 and the 
missile 10. 

[0050] The coordinate aides module 414 receives as aides the stored navigation 

message and pseudorange derived from the receiving station GPS data in accordance with 
tracking aides, Doppler aides or sky search, and coordinates the received aides with the input 
translator GPS data and outputs a navigation message to the PN Code Generator 422 and the 
PN Code Epoch position to the filter and search module 426. 

[0051] The filter and search module 426 controls the input control module 434 to 

locate flie PN Code Epoch position and to sequence to the proper translator GPS data mput 
samples, for example, at time % so that the local carrier generator provides the signal pair 
I*(ti),Q*(ti) and the translator GPS data input sample I(ti),Q(ti) as inputs. The PN Code 
Generator 422 provides tibe PN code in accordance with the navigation message provided by 
the coordinate aides module 414. 

[0052] The I(ti)Q(ti) input signal, I*(ti)Q*(ti) local carrier signal, and PN code are 

processed by the processing stages 410(1) - 410(m), m order to downconvert and wipe off the 
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C/A code of the translator GPS sample input. The outputs of the processing stages 410(1) - 
410(ni) are provided to the filter and search module 426, which filters the results, searches for 
a peak (or peaks) and determines if the results are satisfactory. If the results are not 
satisfactory, &e filter and search module 426 outputs control signals to at least one of the 
local carrier generator 418, PN code generator 422, stage control module 430 and input 
control module 434. The filter and search module 426 thus may control incrementing of ti for 
sequencing to a next input value, changing the aligimient of the PN code or sequencing to a 
different PN code generated by flie PN code generator 422, or changing the number m for 
increasing or decreasing the number of processing stages 410(s) being used. 
[0053] For example, if the results are extremely poor, the filter and search module 426 

may choose to sequence to a next input data sample set, thus throwing away ttie current set of 
samples (P*N samples); change the PN code alignment; and/or increase m. If the results, 
however, are relatively close to being satisfactory, the filter and search module 426 may 
choose to continue processing the current input data sample, realign the PN code and/or 
decrease m for the next sample set. Typically, as processing proceeds, the results of the filter 
and search 226 improve and m is decreased, reducing processing and thus the processing 
time. 

[0054] FIGS. 6A - C shows experimental results obtained using an exemplary 

receiving station 22 in accordance with the present invention. The experimOTt was performed 
with At = 10 Kisec; N=512; chiip=0; Fd=3125.0 HZ; A^195.3125 HZ; and T=5.12 msec/FFT 
window. FIG. 6 A shows a plot of the in-phase signal and noise versus time. FIG. 6B shows a 
plot of the quadrature signal and noise versus time. FIG. 6C shows a plot of spectrum of 
signal and noise versus fi-equency, where C/NO=34 dB-Hz. One can see in FIG. 6C a clearly 
distinguishable peak, providing information that is not readily apparent in FIGS. 6A and B. 
10055] In summary, analysis of GPS data obtained Scorn a fixed receiving station 

approximately simultaneous with GPS data obtained firom a translator mounted within a 
missile is performed using Fourier Transform processing in a series of processing stages. 
Continual tracking of the data is performed by making adjustments to phase and rate control 
of the PN code in the receiving station GPS data while reducing the number of processing 



15 



wo 02/50561 PCT/USOl/47908 



Stages. High precision processing resvdts are possible while minimizing processing steps and 
processing time. Additionally, since the processing involves mathematical computations, the 
entire GPS receiver tracking system can be implemented using only software. Furthermore, 
the ability to continually track while making adjustments in accordance with the control 
signals transmitted by the filter and search module 426 is useful for trackmg in conditions 
such as multi-path GPS signals and reflected GPS signals. 

[00561 What has been described herein is merely illustrative of the application of the 

principles of the present invention. Other arrangements and methods may be implemented by 
those skilled in the art without departing ftom flie scope and spirit of the invention. 
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WHAT IS CLAIMED 15;: 

1 . A receiving station tracking system for receiving GPS data and tracking a 
mobile vehicle, the receiving station tracking system comprising: 

at least one processing stage receiving and processing at least a portion of the GPS 
data; and 

a control module receiving and filtering the processed GPS data and ou^utting at least 
one control signal in accordance with filtering results, wherein the at least one control signal 
controls decreasing or increasing a quantity of processing stages of the at least one processing 
stage. 

2. The receiving station tracking system of Claim 1 , wherein the received GPS 
data includes receiving station GPS data received by a receiver associated with the receiving 
station tracking system and vehicle GPS data received by a receiver associated with the 
vehicle, where the receiving station GPS data and the vehicle GPS data are received 
approximately simultaneously by the receiving station tracking system and flie vehicle, 
respectively. 

3 . The recei vmg station tracking system of Claim 1 , wherem the receiving 
station tracking system processes the received GPS data and recoids the processed GPS data 
as time ordered receiving station GPS data and vehicle GPS data samples. 

4. The receiving station tracking system, of Claim 3, wherein the at least one 
control signal controls sequencing to a subsequent vehicle GPS data sample of fte GPS data 
samples. 
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1 5. The receiving station tracking system of Claim 1 , wherein a respective 

2 processing stage of the at least one processing stage includes processing means for performing 

3 matched filter processing, fest Fourier Transform processing, and square and sum processing 

4 on the vehicle GPS data; and 

5 wherein an output of flie processing means of the respective processing stage is 

6 provided to the control module. 

1 6. The recdving station tracldng system of Qaim 1 , wherein the receiving station 

2 tracldng system is implemented in software having a plurality of programmable instructions 

3 configured for execution by the at least one processor. 

1 7. The receiving station tracking system of Claim 2, wherdn the receiver 

2 associated with the vehicle is included in a translator mounted on the vehicle, and tiie 

3 translator transmits signals including the vehicle GPS data to the receiving station tracking 

4 system. 

1 8. The receiving station tracking system of Claim 2, wherein the receiving station 

2 tracking system further comprises: 

3 a PN code generator module providing to the at least one processing stage a FN code 

4 based on a navigation message of the receivmg station GPS data; and 

5 wherein the at least one control signal furflier controls the PN code generator module 

6 for controlling at least the phase and rate of the PN code. 

1 9. The receiving station tracking system of Claun 8, wherein each of the 

2 processing stages of the at least one processing stage processes the PN c^de delayed by a 

3 respective delay. 
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10. The receiving station tracking system of Claim 2, further comprising a local 
carrier generator module for generating and providing to the at least one processing stage at 
least one periodic signal; and 

wherein at least one periodic signal associated with the vehicle GPS data is provided 
to the at least one processing stage. 

1 1 . The receiving station tracking system of Claim 1 0, wherein the at least one 
control signal controls the local carrier generator module for controlling at least the rate of the 
at least one periodic signal. 

12. The receiving station tracking system of Claim 2, wherein the GPS data 
includes data extracted from the receiving station GPS data. 

13. The receiving station tracking system of Claim 12, wherein the extracted data 
includes at least navigation message data and pseudorange data. 

14. The receiving station tracking system of Qaim 12, wherein the receiving 
station tracking system further includes a coordination module for providing timing and 
polarity data for coordmating the extracted data with the vehicle GPS data. 

15. The receiving station tracking system of Claim 14, wherein the receiving 
station tracking system further includes a peak interpolation module receiving the timing data 
from the coordination module, wherein the peak iiiterpolation module performs a peak 
interpolation algoriflmi on the filtered data, processes the results of the peak interpolation 
algorithm with the timing data, and outputs tracking data associated with the vehicle GPS 
data. 
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1 16. A receiving station tmcldng system for receiving GPS data and tracking a 

2 vehicle, the receiving station comprising: 

3 at least one processing stage receiving and processing at least a portion of the GPS 

4 data; 

5 a PN code generator module providing to liie at least one processing stage a reference 

6 PN code; and 

7 a control module receiving and filtering the processed GPS data and outputting at least one 

8 control signal in accordance with filtering results, wherein the at least one control signal 

9 controls the phase and rate of tiie reference PN code. 

1 17. The receiving station tracking system of Qaim 16, wherein the received GPS 

2 data includes receiving station GPS data recdved by a receiver associated with the receiving 

3 station and vehicle GPS data received by a receiver associated with the vehicle, wherein the 

4 receiving station GPS data and the vehicle GPS data are received approximately 

5 simultaneously by the receiving station tracking system and vehicle, respectively; and 

6 wherein the reference PN code is derived &om the receiving station GPS data. 

1 18. The receiving station tracking system of Claim 16, wherein each of the 

2 processing stages of the at least one processing stage processes the reference PN code delayed 

3 by a different delay. 

1 19. The receiving station tracking system of Qaim 16, wherein the at least one 

2 control signal controls decreasing or increasing a quantity of processing stages of the at least 

3 one processing stage. 
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1 20, The receiving station tracking system of Claim 17, wherein the receiving 

2 station tracking system processes the received GPS data and records the processed GPS data 

3 as time ordered receiving station GPS data and vehicle GPS data samples. 

1 21 , The receiving station tracking system of Claim 17, wherein a PN code 

2 associated with the receiving station GPS data is based on a navigation message of the 

3 receiving station GPS data, 

1 22. A method for tracking a vehicle by a receiving station receiving GPS data 

2 comprising the steps of: 

3 processing at least a portion of the received GPS data by a series of 

4 processing stages; 

5 filtering the processed GPS data to obtain filtered GPS data; 

6 controlling subsequent processing of the processed GPS data based on the filtered GPS data, 

7 wherein controlling includes controlling an amount of processing stages of the series of 

8 processing stages to be subsequently used for processing the received GPS data, 

1 23 . The mefliod of Claim 22, wherein the GPS data includes vehicle GPS data 

2 received by the vehicle and provided to the receiving station and receiving station GPS data 

3 received by the receiving station, wherein the vehicle and receiving station GPS data are 

4 received approximately simultaneousljr, 

5 wherein the step of continually processing fiirther includes the step of processing tiie 

6 vehicle GPS data in conjxmction with a PN code of the receiving station GPS data; and 

7 wherein the step of controlling fiirflier comprises the step of controlling at least the 

8 phase and rate of the PN code. 
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